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ABSTRACT 


An approach to large-eddy simulation (LES) is developed 
whose subgrid-scale model incorporates filtering in the time 
domain, in contrast to conventional approaches, which ex- 
ploit spatial filtering. The method is demonstrated in the 
simulation of a heated, compressible, axisymmetric jet, and 
results are compared with those obtained from fully resolved 
direct numerical simulation. The present approach was, in 
fact, motivated by the jet-flow problem and the desire to ma- 
nipulate the flow by localized (point) sources for the purposes 
of noise suppression. Time-domain filtering appears to be 
more consistent with the modeling of point sources; more- 
over, time-domain filtering may resolve some fundamental ih- 
consistencies associated with conventional space-filtered LES 
approadies. 


1 INTRODUCTION 

By definition, direct numerical simulation (DNS) is the nu- 
merical solution of the Navier-Stokes equations without re- 
course to empirical models. In concept, the fluid motions are 
resolved down to the Kolmogorov length scale, at which ed- 
dies succumb to viscous dissipation. Consequently, for high 
Reynolds number flow, the computational requirements of 
fully resolved DNS are staggering. 

In contrast, in large-eddy simulation (LES), the large scales 
of motion are resolved in space and time on a suitable compu- 
tational grid; however, the effects of the subgrid-scale motions 
on the evolution of the large scales are modeled. Relative to 
DNS, LES is conducted on relatively coarse grids at reason- 
able computational expense. In practice, LES involves filter- 
ing the'Navier-Stokes equations in space or time or both. The 
filtered equations of motion contain subgrid-scale (residual) 
stress terms whose effects must be modeled. 

Both DNS and LES can be classified fundamentally as tem- 
porally or spatially evolving. The distinction between tern- 
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poral and spatial approaches is muddied by most applica- 
tions of LES. In particular, nearly all current LES approaches, 
whether temporal or spatial, exploit spatial filtering. In their 
review paper, Moin and Jimenez (1993) state: “Historically, 
temporal filtering has not been used," one presumes for rea- 
sons of computational efficiency. However, because the roles 
of space and time are fundamentally interchanged in tempo- 
ral vs. spatial simulations, we suggest that spatial filtering 
is more appropriate for temporal LES, and conversely, tem- 
poral filtering is more appropriate for spatial LES. Indeed, 
time-domain filtering may remove some of the conceptual and 
practical inconsistencies that have been observed by practi- 
tioners of LES, a few of which are discussed briefly below. 

First, as intimated by Germaoo (1992), LES can be viewed 
as lying Somewhere near the middle of a spectrum of numer- 
ical solution techniques with DNS at one end and Reynolds- 
averaged Navier-Stokes (RANS) at the other end. In our opin- 
ion, this point of view is most self-consistent if thpe-domain 
filtering is exploited in LES as it is in RANS. Secbnd, Moin 
and Jimenez (1993) observe that the operations of filtering 
and differentiation do not eonimirte on a non-uniform mesh. 
Consequently, most sub^rid-scale models inadvertently im- 
pose different levels of dissipation in different regions of the 
computational domain, a problem made worse on the highly 
stretched grids associated with complicated geometries. This 
problem should be circumvented by temporal filtering in con- 
junction with uniform time increments. Third, again accord- 
ing to Moin and Jimenez (1993): “In LES, it is highly de- 
sirable for the filter width to be significantly larger than the 
computational mesh to separate the numerical and model- 
ing errors. Practical considerations, however, usually require 
the filter width and mesh to be of the same order. In this 
case, there does not appear to be a necessity for higher than 
second order numerical methods for LES.” In contrast, for 
the present temporally filtered approach, the filter width is 
typically an order of magnitude larger than the time step. 
Fourth, it may be desirable in spatial DNS or LES of certain 
physical problems (e.g., jet flow; to allow for time-dependent 
localized (point) sources as a means of manipulating the flow 
for the purposes of control. For example, such sources could 
be used to introduce local disturbances to enhance or inhibit 
mixing. Dakhoul and Bedford (1986) suggest that spatial fil- 
tering is fundamentally inconsistent with the introduction of 
point sources, whereas temporal filtering of a point source is 
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well defined. 

Whereas Dakhoul and Bedford (198fi),aii4^1dama (1990) 
propose and develdp space-time filters forL®S; the author is 
unaware of any purely time-filtered approach. In the next few 
sections, we develop and demonstrate a spatial LES concept 
based on filtering in the time domain, and we apply the ap- 
proach to the investigation of large coherent structures (CS) 
in a heated subsonic axisymmetric jet. For several reasons, 
the jet-flow problem is well suited to the particular LES ap- 
proach. First, free shear layers, jets, and wakes, whose mean 
streamwise velocity profiles are inflectional* are myisddly un- 
stable to disturbances of a broad spectrum of frequencies 
(wavelengths). As a consequence, DNS of three-dimensional 
(3D) unbounded shear flows is presently impractical because 
of tne extremely fine grid resolution required, and some sort 
of subgrid-scale dissipation is virtually a necessity. Second, 
the problem is of immediate practical interest to the field of 
computational aeroacoustics (CAA). Specifically, it Is gener- 
ally believed that, for supersonic jets, most of the noise origi- 
nates from the CS rather than from the small-scale turbulence 
(Seiner, 1984, and Tam, 1995). For subsonic jets, the origin 
of mm i*less certain. In his recent review paper, T^Pi (1995) 
expresaei tbe view that subsonic jet noiifpri^nffoe primar- 
ily from fine-scale turbulence. Our results fpwliminary at the 
present time) would suggest otherwise. Thus, spatial LES 
may provide a tool by which to investigate the physics of 
noise production and suppression ia jets. Last, because no 
watts are present, we avoid for the time being the difficulties 
experienced by many subgrid-scale models in the vicinity of 
sohd boundaries. 


The input to Eq. 1 is the raw signal s(t), and the output is 
the continuous filtered signal, denoted by 3(t, A). From ele- 
mentary calculus, the followup property of the filter defined 
by Eq. 1 is readily derived: 

s{t, 0) = Urn s(t , A) = s(t) (2) 


In general, I / 3 If A represents some moderately large tem- 
poral window, then filtering «(£) via Eq. 1 wUl tend to remove 
oscillations of high frequency relative to A while preserving 
low-frequency oscillations, which defines a “low-pass” time- 
domain filter. 


MNgARB, 


AUSAL FILTERS: Let us digitize 


the continuous signal $(f j such that s» = «(**), where U = iAt, 
and At is the (constant) time interval between samples. Typ- 
ically, for applications to LES, A should be an order of mag- 
nitude larger th an At. The approximation of Eq. 1 by a linear 
quadrature rule results in its discrete analog 


m 

*» = yipj^i-j 

i=0 


(3) 


where the filter coefficients pj are determined to give the fil- 
ter certain desirable properties (e.g., low-pass characteristics, 
stability, and high-order accuracy at low frequencies) . Follow- 
ing Press et si. (1986), we generalize the linear digital filter 
given in Eq. 3 to allow the use of previously filtered data. 
Specifically, suppose 


In the next' action, time-domain filtering is discussed in 
general, and a prototype causal digital filter is developed. The 
governing equations for DNS are discussed in fection 3, and 
the equations Sire modified for LES based on adAptation of the 
so-caued SEZHu (“says who”) model ofSpezlaleand cowork- 
ers (1988). The present model differs from the $E3&u model 
in that i| exploits temporal rather than aptfe! fitef e jjwg The 
axisymmetric-jet test case is defined iq j@i£ibn A:Tne nu- 
merical approach to the solution of the governing equations, 
adapted from Pruett et al. (1995), is addres se d briefly in Sec- 
tion jjht Results of LES for the axisymmetric-jet problem are 
presentadin Section 6 and are compared with weU-reeolved 
DNS results. Finally, some brief conclusions are offered in 
Section 7. 


2 CAUSAL FILTERING 


Time-domain fitters foil into either of two categories: causal 
or acausaL Fbr application to LES, only causal filtering is re- 
alisable* By definition, causal temporal filters exploit present 
and pastdainrination only, the future being inaccessible. Con- 
sequently, In this section, we consider prototypical c ontinu ous 
and discrete causal filters. A continuous filter is presented for 
conceptual purposes; an analogous discrete filter is exploited 


in practice. 


ij* if 


IUER: If sft) represents a 
Cptitinubug signal in rime then a COntinUOUB lOW- 
pass causal filter can be constructed simply by integrating 
the signal over the interval A, the temporal window width, as 
follows: 


,(, ' a) =s£ 


s(r)dr 


( 1 ) 


= 2 P ’ 8i - j + S qkli ~ k ( 4 ) 

$=> 0 Mml 

The filter of Eq. 4 is “nonrecursive” if 7 * = 0 for all k and 
“recursive" if, for at least one k, g* ^ 0, in which case the 
current value of the filtered quantity is a linear combination 
of previous unfiltered and filtered values. 

FREQUENCY RESPONSE: It is instructive to examine 
the frequency response of the filter associated' bfith Eq. 4. 
From Press et al. (1986), the transfer function, which quan- 
tifies the frequency response, is given by 


I - E*»i 


( 5 ) 


where t = >/-T ft = w*Af* is the dimensionless frequency, 
w* = 2 nf m is the dimensional circular frequency, and /* is the 
dimensional physical frequency. (Throughout this work, we 
denote dimensional quantities by asterisks.) In general, the 
frequency response of a recursive filter is related to a ratio- 
nal polynomial function in the complex variable 1/C, where 
C = e*°. Thus, recursive filters are to nonrecursive filters what 
compact-difference operators are to standard finite-difference 
operators. The rational polynomial form of the transfer func- 
tion allows considerable latitude in shaping the frequency re- 
sponse. Fig. 7 compares the modulus of the transfer ftmctkm 
of a prototypical low-pass digital recursive filter with that of 
an idealized; “spectral cutoff” filter, fbr a nominal cutoff fre- 
quency ft'. Note that, for the spectral cutoff filter, ? = J, 
which, as we have mentioned previously, is not true in gen- 
eral. The transfer function of the digital filter can be made 
to more closely approximate the spectral ideal at the expense 
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of including more and more history (i.e., by using larger and 
larger values of m and n.) 

. . For the purposes of time-filtered LES, the design con- 
straints^^ the discrete filter are: 1) stability for all U\ 2) 
H(0) = 4 ; 3 ) high-order accuracy; 4) |i/(fi)| = 0 for fl > H c ; 
turn 5) as little storage required as possible. For reasons 
to be addressed fully m a subsequent paper, one is lead to 
the fortunate if surprising conclusion that second-order fil- 
ters are optimal for the present application to LES. First, 
second-order causal filters require relatively little storage for 
hiitpiy.’ 'Second, to avoid mixing the truncation errors of the 
filter am numerical method, one should invoke a filter with 
nohighfer order than thst of the time-advancement scheme. 
Thus, a second-order filter is compatible with the present 
third-order time advancement scheme (see Section 5). Third, 
and a stibtle point, it can be shown that second-order filtering 
isoonsiitent with the underlying subgrid-scale model (Pruett, 
1996b) . Consequently, for our purposes, we have followed the 
design procedures outlined in the digital signal filtering text 
by Strum and Kirk (1988). Formally, our prototype filter is a 
second-order, impulse-invariant, digital Butterworth low-pass 
filter, for which m = h — 2, and po — 0. Technically, But- 
terworth filters are “all-pole” filters, whose transfer functions 
are maximally flat in the vicinity of the origin. Unfortunately, 
however, the transfer functions of Butterworth filters do not 
vanish identically for large values of fl. In practice, this is 
not a problem as will be shown subsequently. The nominal 
cutoff frequency for our prototype filter is fl' = 1.0, for which 
|tf(fl' c )| 2 =0.5. 

A suitable generalization from the prototype filter to a 
tunable-cutoff low-pass filter is made by incorporating a pa- 
rameter Re, defined as the ratio of the actual and prototypical 
cutoff frequencies, namely 


_ ftc _ At 
~ fij. — A 


( 6 ) 


Note that Re -*■ 0 for the discrete filter is analogous to A -> oo 
for the continuous filter of Eq. 1. Conversely, as Re -4 oo, 
A -> 0, in which case I(t, A) -> a(t) by the property of Eq. 2. 
Fig. 7 shows the frequency response of the present filter far 
Re = 425, a value typical for the current time-filtered LES 
approach. Note that high-frequency oscillations a?e virtually 
eliminated by the present filter, as desired. 


3 GOVERNING EQUATIONS 

We first specify the governing equations for DNS and then 
present the governing system as modified for LES. 

DNS: As a basis on which to evaluate various LES solutions, 
we require a well-resolved DNS solution for the axisymmetric- 
jet problem. For a compressible fluid, it is appropriate to 
define a, fluid state vector [p,p,T,u,v,w] T comprised of the 
density p, pressure p, temperature T and velocity components 
u, v, and to. The governing equations for the axisymmetric- 
jet problem are adapted from those presented in Pruett et al. 
(1995) for a body-fitted coordinate system x = |x, 0, z] T on 
an axisymmetric body, where z as the arc length along the 
body, 6 is the azimuthal angle, z is the coordinate normal to 
the body, r = R+zcos^is the radial coordinate, R(x) is the 
body radius, anid 4>(x) is the angle of the surface tangent to 
the body. For the jet-flow application, R = </> = 0, in which 


case the coordinate system degenerates to : = r with x as the 
axial coordinate, .whereby u and w become the axial and ra- 
dial velocities, ~ ^fesSpectively. Because R = 0, the equations are 
geometrically singular along the jet axis. For the continuity 
equation, the singularity is removed by applying L’Hopital’s 
rule along the axis (z = 0). The singularity is not problematic 
for the momentum and energy equations because axial bound- 
ary conditions replace the governing equations along the axis. 
Specifically, considerations of symmetry require that the az- 
imuthal velocity (v) vanish everywhere and that 


dT 

dz 



= 0 


(* = 0 ) 


( 7 ) 


LES AND SUBGRID-SCALE MODEL: If the compress- 
ible Navier-Stokes equations (CNSE) are filtered in the time 
domain according to Eq. 1, the resulting equation system is 
formally identical to that of Eqs. (15), (16), and (34) of Er- 
lebacher et al. (1992), where overbars and tildes distinguish 
conventionally filtered and Favre-filtered quantities, respec- 
tively. In general, the use of Favre-filtered (density-weighted) 
variables reduces the complexity of the filtered CNSE. Specif- 
ically, for example, the Favre-filtered axial velocity is defined 
as 

u = ^ (or pu = pit) (8) 


Other Favre-filtered quantities are defined analogously. For 
the filtered equations, the fluid state vector is comprised 
of a mixture of conventionally and Favre-filtered quantities, 
namely [p,p,T, u, ©.iS] 1 *. The filtered governing equations 
contain residual stresses not present in the origins! equations, 
which are decomposed into Leon&rd-Stress, cross-stress, and 
Reynolds-stress terms denoted by L, C, and R, respectively, 
following the notation of Erlebacher et al. (1992). Of these, 
G and R must be modeled; L can be computed. 

Several candidate subgrid-scale models are available; how- 
ever, for our present purpose, we adapt the SEZHu model 
(Speziale et. al, 1988) as adapted by Erlebacher et al. (1992). 
Several considerations favor this selection. First, to demon- 
strate the time-filtered approach, the filtering process must 
participate in the model, rather than simply serving as a con- 
ceptual framework as it does In some models (e.g., Smagorin- 
sky). Second, although present results suggest that a dy- 
namic model (Germano et al., 1991) is desirable, we wanted 
initially to avoid some of the pitfalls of dynamic subgrid-scale 
models, particularly, the need to smooth the computed model 
constants. Finally, the SEZHu model Is extremely well docu- 
mented in Erlebacher et al. (1992). 


Formally, our implementation of the SEZHu model is vir- 
tually identical to that of Erlebacher et al. (1992), except, 
of course, that the filter is temporal; In tensor notation, the 
dimensionless filtered equations, with the modeled terms de- 
noted by underlines, are 

yM 2 p s pf (9) 


ap g(gufc) _ 

at dx k 


( 10 ) 


at 


+ [/JUjfeUj + p(ufcfij - dfctq)] = 


ap_ 

dx k 


+ ^[ ( A**+ff t)S U ] 
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dp. JL 

dt dxk 


-(7-1 )PD + 


where 


put + - iikf) 


d 

dx k 






D = 


dui, 

Ox* 


is the resolved-scale dilatation, 


( 12 ) 


(13) 


S u = 2 (e*i - - DSu ) 


(14) 


Su is the kronecker delta, and Su is the resolved-scale strain- 
rate tensor, namely 


e« = 2 


, dvi 

dxi dxk J 


(15) 


For brevity, the physical viscosity and thermal conductivity 
are denoted, respectively, as 


» = JL 

^ Re 


= 


M 2 RePr 


(16) 


where Re, Pr, and M axe the dimensionless Reynolds, 
Prandtl, and Mach numbers, respectively. Similarly, the eddy 
viscosity and the eddy thermal conductivity are given by 


HX m CrPjfll 1 / 2 


k t = 


fi T 


•yM a Prr 


(17) 


where l is a length scale to be defined shortly, Pry is. the 
turbulent Prandtl number, 7 is the ratio of specific heats, 

II = SuSu (18) 

Whereas the underlined terms on the right-hand sides of the 
governing equations model R, the underlined terms on the 
left-hand sides, which are properly referred to as the resolved 
stresses, are computable by filtering the resolved fields. For 
the present subgrid-scale model, the resolved stresses model 
the sum L + C. 

In the original SEZHu subgrid-scale ihodel (Spesdate et al., 
1988), the Reynolds stresses were split into deviatoric and 
isotropic parts, which ware modeled separately. In the more 
detailed paper of Erlebacher et al. (1992), on which the 
present approach is based, the isotropic part is dis reg a r ded 
on the rationale that its contribution should be small ror tur- 
bulent Mach numbers M t < 0.6 (a constraint satisfied by most 
compressible flows). Several other subtleties of the ir _ 
tation of the model are not immediately apparent upon the 
study of Erlebacher et al. (1992). First, the viscous-stress 
terms of the filtered momentum equations, and the dissipa- 
tion function 8 and thermal-stress terms of the filtered en- 
ergy equation, formally involve conventionally filtered rather 
tha n Favre-filtered quantities. Because these quantities are 
unavailable, however, they are approximated by their Favre- 
filtered equivalents. Second, terms that arise from subgrid- 
scale fluctuations of fh> and k v are neglected. Third, the re- 
solved stresses are computed using p rather than p, the latter 
of which is unavailable. These approximations should be con- 
sidered as additional modeling errors. 


MODEL CONSTANTS: Equation 17 requires values for 
three constants. Following Erlebacher et al. (1992), we use 
Prr = 0.5 and. C T = 0 . 012 . It remains to determine l, which, 
for the original SEZHu model, is a characteristic length scale 
related to spatial grid resolution. Specifically, Erlebacher et 
al. (1992) show that If = cAx is optimal (in the sense of pre- 
serving the Galilean invariance of certain terms of the model) 
for c * 2 , where Ax is the afctual computational grid spacing. 

Here, we must determine l based on the choice of the tempo- 
ral scale A (or equivalently, Re), for which purposes we appeal 
to results mean toe area of hydrodynamic stability. From lin- 
ear stability theory, we know that jets, wakes, and free shear 
layers are dispersive; i.e., wavesofdifferent frequencies prop- 
agate Jit different phase velocities. However, disturbances of 
moderate to high frequencies propagate at a velocity approx- 
imately that of Dw, the average of the jet and ambient ve- 
locities. Accordingly, we define Az* as the characteristic size 
of an eddy associated with a disturbance of cutoff frequency 
/*; that is 


a _• _ U& v 2*UgyAt* 

S: Rc 


(19) 


Finally, we note that the CNSE are recovered from the 
governing equations in the limit as A -4 0. 


4 TEST CASE 


The numerical test case was chosen to approximately replicate 
an acoustics experiment that is being conducted at NASA 
Langley Research Center. Specifically, we investigate a heated 
subsonic (At m, 0 . 8 ) jet exhausting into a nearly quiescent at- 
mosphere- The jet temperature Tj , on which the Mato num- 
ber is baaed} is 60Q F (10S9.6 R), and toe ambient tempera- 
ture is T0 F (529.6 R). The nominal jet radius is Rj = 0.5 in. 
(0.0417 ft j. The ambient prestore in the physical experiment 
is approximately one atmosphere (2HK) ps£). However, this 
results in too high a Reynolds numher for a DNS computa- 
tion of reasonable expense; consequently, the computational 
experiment assumes an ambient pressure 10 percent that of 
the physical experiment; that is, 216 psf, Whkh results hi 
Re s 10153 based on the jet conditions and the nominal jet 
radius. In the physical experiment, the ambient air was quies- 
cent. However, computational experiments with unbounded 
shear layers typically encounter numerical difficulties (as did 
the present work) whenever the a mbient stream is perfectly 
qu i esce n t (Tannehili et al., 1984). Consequently, it is cus- 
toms^ for the jet to exhaust into a coflowing stream with a 
velocity of a few percent of the jet velocity. For the present 
problem, we use an ambient to jet velocity (IT?) ratio of 10 
percent. We further assume that the jet is folly expanded, in 
which case, in the absence of any disturbances, the pressure 
is constant both radially and axially. 

la toe governing equations and in toe results to follow, all 
lengths hate been normalized by Jtjj and toe velocities, tem- 
perature, and density, have been normalized by Iff, Tf, and 
Pj, respectively. Pressure is normalized by P)Uf 2 . 
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5 NUMERICAL METHODOLOGY 

Spatial DNS and LES can be viewed as three-step processes. 
First, an unperturbed time- independent base state is ob- 
tained, usually by boundary-layer techniques. Second, the 
base state is Subjected to either random or temporally peri- 
odic perturbations, which are typically imposed at or near the 
computational inflow boundary. The structure of these dis- 
turbances is commonly obtained from linear stability theory, 
or mote recently, from parabolized stability equation (PSE) 
methodology. Third, the Spatial evolution of the propagating 
-disturbances is computed by solution of the complete Navier- 
Stokes equations, with (LES) or without (DNS) subgrid-scale 
models. 1 ' We discuss in turn each of these steps in the context 
of the Current problem. 

• COMPUTATION OF THE BASE STATE: The applies 
tion oL standard fully implicit boundary-layer techniques to 
the agUymmetric jet revealed an unanticipated computational 
difficult; namely, the Jacobian matrix associated with the 
iteration procedure was extremely ill-conditioned and the it- 
eration} did not converge. The computational problem arises 
from a -reversal of sign in the transverse velocity experienced 
■by internal boundary-layer flows such as jets, wakes, and free 
shear layers. To circumvent this numerical difficulty, a semi- 
implicit boundary-layer technique was developed, which in is 
documented in Pruett (1996a). The interested reader is re- 
ferred to this paper for details. Also to avoid numerical diffi- 
. cultie$£the internal shear layer is given finite thickness at the 
lip of the jet. The boundary-fryer solution can be viewed as 
an unstable equilibrium state of the CNSE. 

, 5ITION OF THE DISTURBANCES: The nature of 
.... .. njuties is different for wall-bounded and free-shear flows. 
Specifically, wall-bounded flows are subject to viscous insta- 
bilities, for which typically only a relatively narrow band of 
frequencies are unstable. In contrast, free-shear layers, jets, 
and wakes are subject to inviscid instabilities over a broad 
range of frequencies. Relative to viscous instabilities, inviscid 
instabilities experience rapid growth rates. Thus, for simulat- 
ing instability waves in wall-bounded flows, it is essential that 
the imposed disturbances be consistent with eigenfunctions 
obtained from stability theory; otherwise, one introduces spa- 
tial transients that may corrupt the particular instability of 
interest. On the other hand, one can be somewhat cavalier in 
imposing disturbances in free-shear flows because of the flow’s 
tendency to rapidly organize arbitrary disturbances into the 
dominant eigenmodes. Consequently, following Mankbadi et 
al. (19941, at the inflow boundary, we impose a temporally 
periodic fluctuation comprised of a few harmonics of specified 
frequencies, but whose structure is not derived from stability 
theory. At present, we impose the disturbance only through 
the stre&mwiae velocity. Specifically, at the inflow boundary 
X » 2?0 


u(t,x Q ,z) 

= ub(x 0 ,z) + eu'(t, z) 

(20) 


— 4>(z) [sin(u;/f) + cos(0.5u;/t)] 

(21) 

m 

= exp{— (2(s — l)) 4 ] 

(22) 


where the subscript B denotes the base state, the prime de- 
notes a fluctuating quantity, and the subscript / denotes the 
fundamental frequency* The function <f>(z) is used to shape 
the disturbance profile so that the disturbance is largest near 
the edge of the jet but essentially vanishes along the jet axis 
and at the far-field boundary. Numerical experimentation re- 
veals the most rapid development of the jet for u/ = tt, which 


corresponds to a Strouhal number (St = f)Rj/Uj) of 0.5, 
in keeping with the observations of Mankbadi et al. (1994). 
Following the early computational investigation of free-shear 
layers by Riley afid Metcalfe (1980), we include an out-of- 
phase subharmonic component to enhance pairingof adjacent 
vortices. At present, we use a forcing amplitude of e « 0.005, 
which is ramped up slowly (over a time interval of one pe- 
riod of oscillation at the fundamental frequency) to minimize 
temporal transients. 

DNS AND LES METHODOLOGIES: For both the DNS 
and LES, we adapt the high-order numerical scheme of Pruett 
et al. (1995), to which the reader is referred for .details. 
Briefly, this algorithm exploits fully explicit time advance- 
ment, high-order compact-difference methods (Lele, 1992) 
for aperiodic spatial dimensions, and spectral collocation 
methods for periodic spatial dimensions. Spedfitklly, for 
the present axisymmetric-jet application, we use fourth- and 
sixth-order compact difference schemes in the axial and ra- 
dial dimensions, respectively. The azimuthal dimension, of 
course, does not come into play for the axisymmetric case. 
The method of Pruett et al. (1995) uses a variable step for 
time advancement in the context of a three-stage, low-storage 
Runge-Kutta (RK3) scheme. However, the present LES appli- 
cation, which involves temporal filtering, requires a constant 
time step. Consequently, the original Runge-Kutta tempo- 
ral integration has been replaced by a fixed-length, multiple- 
step, third-order Adams-Bashforth (AB3) technique. An ad- 
ditional motivation for replacing the RK3 method was that it 
was not immediately clear to the author how temporal filter- 
ing would interact with time advancement whenever multiple 
stages per time step were involved, the fear being the pos- 
sibility of numerical instability. The storage requirement for 
the algorithm with AB3 time advancement is about 150 per- 
cent that of the original algorithm with RK3. In general, 
multiple-step methods are not self starting. The AB3 inte- 
gration is started initially with a single first-order Eider step 
followed by one second-order Adams-Bashforth step. Because 
the perturbation is ramped slowly to full amplitude, and the 
Initial state is in (near) equilibrium, the initial loss of tempo- 
ral accuracy is inconsequential. 

For both the DNS and the LES, the symmetry conditions 
riven by Eq. 7 are imposed along the jet axis. At the in- 
flow boundary, for the present; arisymmetric-jet- problem, the 
Jbw is everywhere subsonic, and one characteristic points up- 
stream. Consequently, not all flow variables can be specified. 
Currently, we specify v, w, T, and the incoming Riemann in- 
variants. At the far-field boundary (z — zmax), we adapt 
tiie irifi-teflecting boundary conditions of Thompson (1987) 
as Modified by Pruett et al. (1995). At the outflow boundary, 
weeriptoit a buffer-domain approach (Streett and Macaraeg, 
1989/1990). Near the outflbw boundary, a buffer zohe of fi- 
nite width is constructed in which both the base state and 
the governing equations are modified to ensure that all waves 
propagate out of the domain. 

Finally, we note that the present LES algorithm is one of 
few to incorporate high-order numerical methods, another be- 
ing that of El-Hady and coworkers [5]. 

COMPUTATIONAL EFFICIENCY: On the same grid, an 
LES computation with the present algorithm requires not 
qwte twice the computational effort as DNS and approxi- 
mately 2.5 times the storage for 3D flows (twice the stor- 
age for two-dimensional or axisymmetric flows). Most of the 
additional memory is relegated to storage of the time histo- 
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ties of quantities associated with the time-filteredaj>proach. 
One must keep in mind, however, that, by definition; LES 
allows computations oncoarser grids than DNS. If; for the 
present LES algorithm, fei example, the grid reablutloB rela- 
tive to DNS could be reduced by a factor of threeineach of 
the three spatial dimensions and time, then storage require- 
ments would be diminished by a factor of appeoxknstely ten, 
and processor time would diminish by a factor of approxi- 
mately 40. Thus, measures of efficiency in LES rausteofisider 
not only nominal storage aqd operation counts, but also the 
potential grid-coarseaifig fector, which could conceivably be 
higher for temporally filtered LES than for conventional ap- 
proaches. ' 


6 RESULTS 


For the DNS and LES results presented below, the computar 
tional domain was 


0 < * < 20 


0 < z < 5 


1 v 


(23) 


The length of the domain was sufficient to allow one pairing of 
the adjacent vortices shed at the edge of the jet^Tke final 16 


main; results within the buffer domain should be disregarded 
as unphysical. For convenience in presmiting results, we de- 
fine tp, the time in periods of oscillation at the fundamental 
disturbance frequency. 

DNS: The DNSresults were obtained on the computational 
domain defined above at an extremely fine spatial grid reso- 
lution of 1280 x 512 and a temporal resolution of 2048 steps 
per (fundamental) disturbance period. To arrive atthis res- 
olution, computations were made on successivehr finer grids 
beginning from a coarse. grid of 256 x 128. For each spatial res- 
olution, an estimate bf the temporal resolution necessary for 
stability was made based on stability analyses of model aavec- 
tion and diffusion equations. All computations except fe* that 
on the finest grid eventually “blew up” due to numericafinsta- 
bilities associated with unresolved scales. In contrast, Fig. 7 
shows instantaneous contours of constant density at tp = 18 
for the fully resolved computation. On the finest grid, the 
DNS computation required in excess of 20 CPU hours on a 
■Cray C90. Consequently, calculations at a higher Reynolds 
number would have been impractical given the computational 
resources available. 

The present fine-grid DNS represents one of the most accu- 
rate and highly refined computations of an unbounded shear 
layer of which we are aware, another being that of Coburns et 
al. (1995), bod Fig. 7 affords considerable detail. If is inter- 
esting to note a striking correlation between the contours of 
constant density and those of constant vortkity (not shown 
due to space limitations). Both quantities clearlyshow the 
roll-up of the shear layer at the jet's edge into a vortex street 
and the subsequent pairing of adjacent vortices, phenomena 
common to unbounded shear flows. A similar comparison of 
vorticity and pressure contours is also most revealfrig. Not un- 
expectedly, the centers of low pressure correspond precisely 
with the centers of the large vortices. As adjacent vortices 
merge, the individual pressure lows are replaced by lareer 
and stronger low pressure regions. High pressure regions lie 
between adjacent vortices. 

T.F.S: For LES of jet flow, the trick is to find an appro- 
priate amount of eddy viscosity. If the LES subgrid-scale 


model is insuffiCksitiy dissipative, the computation 
up. On the other hand, if the model is excessively dissipative, 
the instabilities that result in vortex shedding and pairing 
are suppressed or diminished. It appears possible, however, 
that an intermediate amount of dissipation wifi preserve the 
large-scale. features of the flow while preventing numerical^* 
stabilities associated, with, unresolved scales. In particular, 
Figs, 7 and 7 presentinstantaneous contours of constaatden- 
sity and pressure, respectively, obtained from an LES compu- 
tation wkh ^ spatial grid resolution of 432 x 192, coarser by a 
factor of appnsdmatriy three in each direction than the DNS 
computation presented previously. The time is tp = 18, and 
the temporalresolutjon is 1024 time steps per period. Be- 
cause fully explicit time-advancement schemes typically yield 
over-resolution in time, practicality demands a filter width A 
that is approximately, say, an order of magnitude larger than 
the time step. Here, we used A/At = 8 (ft* = 0.125), At this 
cutoff frequency, the maximum eddy viscosity was approxi- 
mately eight times the maximum physical viscosity. Whereas 
the DNS calculation required 20 CPU hours, the 432 x 192 
LES cakufr&on tequired two CPU hours. (Calculations on 
the coarsest grid ot 256 x 128, for which the computation 
eventually blew up, required only a matter of minutes.) Rel- 
ative to the DNS results of Fig. 7, the shear-layer roll-up 
and pairing events of the moderately dissipative LES com- 
putation are retarded, but ndt prevented. Consequently, we 
believe thfr dKHferately resolved LES could serve as a compu- 
tationalplaiform for the investigation of jet noise. To rids 
end and following Cobnius et al. [3], we extract the com- 
pressible dilatation from the numerical solution; The instan- 
taneous dilatation field at tp = 18 of the moderately resolved 
LES computation is shown in Fig. 7. It would appear that, 
in tezmspftjtedilatation, each Targe vortex appears,?? an 
acoustic quadrupole. These results suggest that significant 
acoustic radiation is associated with the large-scale vortices, 
in contrast to the view expressed by Tam (1995) for subsonic 
jets. ;; 

For the present work in progress, our primary objective has 
been to demonstrate the feasibility and practicality of time- 
filtered spatial LES. As one possible measure of success, we 
examine the resolved subgrid-scale stresses computed by the 
time-filtered subgrid- scale model, the principal component of 
which is pr es ented ha Fig: 7. For an LES computation, rite 
magnitude of the resolved stresses can be viewed as a mea- 
sure of ill resolution, or equivalently, as identifying the re- 
gions where additional dissipation is needed to prevent numer- 
ical instability. Fbr rite present time-filtered approach, these 
terms are well defined, judging by the apparent smoothness of 
the contours in the figure. In summary, these results suggest 
that a time-filtered dynamic subgrid-scale model could be de- 
veloped. We suspect that dynamic modeling (Gerraano et at., 
199!) would improve the present results by limltingtbe eddy 
viscosity only to regions where it is needed. For this reason, 
the dynamic approach is judged to be more appropriate for 
transitional Sows, as was demonstrated for compressible flow, 
for example, in the work of El-Hady et al. (1993). 

In closing, we comment that computations of 2D or ax- 
isynunetric unbounded shear flows are both less and more 
difficult, than computations of 3D flows. Although the total 
computational requirements are lower for simulations in 2D, 
the resolution needed for those two dimensions may well be 
neater than for the same two dimensions of a 3D simulation. 
The reason is that, in 3D, the third dimension provides a path 
for relief of Reynolds stresses that cannot be relieved in 2D. 
Complete validation of the present time-filtered LES concept 
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and its application to aeroacoustics ultimately will require the 
consideration of fully 3D flows. 


7 CONCLUSIONS 

• Tiine-domain Altering for the residual-stress models of 
LBS if a viable concept that should be investigated fur- 
ther. Pre liminar y 'results suggest that time-domain fil- 
tering Stay have significant advantages relative to con- 
vefetioni! space-filtered approaches. 

• The current baseline LES algorithm is one of a very few 
LES algorithms to exploit high-order numerical methods. 

• The present subgrid-scale model, which involves time- 
11 v domain; filtering, might be improved for application to 
i transitional flows by recasting it in a dynamic-model con- 

text (Germano et at., 1991). 

• The present approach to LES appears to be applicable 
to aeroacoustics, as had been hoped. 

• A thorough validation of the time-filtered LES approach 
and its usefulness to aeroacoustics will require LES of 
fully 3D flows. 
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FIGURES 
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Figure 1. Transfer function of prototype second-order 
causal filter compared with spectral-cutoff transfer func- 
tion. 
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Figure 4. Contours of constant pressure obtained from 
LES calculation of resolution 432 x 192 at tp = 18. Con- 
tour levels denote variations of 0.87 < p < 1.26. Buffer 


domain shown. 
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Figure 2. Contours of constant density obtained from 
well-resolved (1280 x 512) DNS calculation at tp = 18. 
Flow is left to right. Jet centerline is along lower bound- 
ary of figure. Contour levels denote variations of 0.82 < 
p < 2.15. Buffer domain shown. 
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Figure 3. Contours of constant density obtained at tp = 
18 from LES of 432 x 192 grid resolution. Maximum eddy 
viscosity is 8 times that of physical viscosity. Relative 
to DNS, vortex roll-up and pairing events are retarded 
but not prevented. Contour levels denote variations of 
0.95 < p < 2.13. Buffer domain shown. 
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Figure 5. Isolevels of dilatation D obtained at tp = 18 
from LES of 432 x 192 grid resolution. Contour levels 
denote variations of —0.10 < D < 0.13. Buffer domain 
not shown. 



Figure 6. Principal component of resolved stresses ob- 
tained at tp = 18 from LES of 432 x 192 grid resolu- 
tion. Contour levels denote variations from -4.2 x 10~ 5 
to 3.3 x 10-*. Buffer domain shown. 
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